function alphaplot3(T,X,radius,volume)
% Plot 3D alpha shape

% Remove inner faces
F = [T(:,2:4); T(:,[1 3 4]); T(:,[1 2 4]); T(:,1:3)];
F = sort(F,2);
F = sortrows(F);
ii = find(~any(diff(F),2));
F([ii;ii+1],:) = [];

% Coordinates
x = X(:,1);
y = X(:,2);
z = X(:,3);

% Plot boundary faces
trisurf(F,x,y,z,'FaceColor','yellow','FaceAlpha',1);
% hold on
% plot3(x,y,z,'k.')
% hold off
axis equal

% Title
str = sprintf('Radius = %g,   Volume = %g',radius,volume);
title(str,'fontsize',12)